IDENTIFICATION OF BOUNDARY CONDITIONS USING 
NATURAL FREQUENCIES IN CASE OF A RING MEMBRANE 
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Abstract 

The problem of finding boundary conditions for fastening of a ring membrane, which are 
inaccessible for direct observation from the natural frequencies of its flexural oscillations, 
O ' is considered. 

Two theorems on the uniqueness of this problem are proved, and a method for estab- 
lishing the unknown conditions for fastening of the membrane to the walls is indicated. 
An approximate formula for determining the unknown conditions is obtained, using first 
three natural frequencies. The method of approximate calculation of unknown boundary 
fN| . conditions, is explained with the help of an example. 

Keywords: Boundary conditions, inverse spectral problem, membrane, natural frequen- 
p '■ cies, Pliicker coordinates, Pliicker relation. 

1. Formulation of the Inverse Problem. 

Let us give the mathematical formulation of the problem, before describing 
the method of its solution. The problem of the radial oscillations of membrane, is 
reduced after making a separation of variables, to the following eigenvalue problem 
: (see [5]) 

y{r) +y{r)/r + X 2 y{r) = 0, (1) 
U 1 (y) = k 1 y(a)-k 2 y(a)=0, U 2 (y) = k 3 y(b) + k 4 y(b) = 0, (2) 

O 

where A, is the spectral parameter, a and b are the short and long radius of the 



ring membrane respectively. 

Now we formulate the inverse of the eigenvalue problem: 



Problem 1 (inverse problem). It is required to find the unknown linear forms 
Ui(y) and U2(y) from the eigenvalues of the problem (l)-(2). 

Let L(y, A) = y(r) + y(r)/r + \ 2 y(r), (U^y), U 2 (y)) T = AX T , 
where A = {aij\i = 1,2; j = 1,2,3,4} is a matrix of order 2x4 such that 
&ii = hi, a\2 = — k 2 , ^23 = a 24 = ^4, &13 = du = a 2 \ = a 22 = 0, and 

X(y) = (y(a),y(a),y(b),y(b)) is a row-vector. Then the eigenvalue problem 
(l)-(2) is equivalent to the following eigenvalue problem 

L(j/,A) = 0, AX T (y) = 0. (3) 

We consider also another eigenvalue problem 

L(y,X)=0 1 AX T (y)=0. (4) 
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In these problems only the matrices of boundary conditions are different. 

Proposition 1. If matrix A is equivalent to matrix A, i.e. there exist non- 
singular matrix S such that A = SA, then problem (4) is equivalent to problem 



The proof is trivial. 

Proposition 2. Matrix A is equivalent to matrix A, if and only if, their 
correspond minors of order 2 are equal, accurate to coefficient. 
The proof is in [3]. 

In the projective geometry we define the Pliicker coordinates of the set of 
equivalent matrices of order m x n to be any one of the equivalent c|^^-tuples 
of its minors of order m. 

Now we shell give the following 

Definition 1. Pliicker coordinates of boundary conditions of the eigenvalue 
problem (3) are called any one of the equivalent 8-tuples of minors (of order 2) of 
its matrix A: (. . . , Aij , . . .). 

From Proposition 1 and Proposition 2 we get 

Proposition 3. // Pliicker coordinates of boundary conditions of the eigen- 
value problems (3) and (4) are identical, then problems (3) and (4) are also 
identical. 

Then Problem 1 is equivalent to the following problem: 

Problem 2. It is required to find Pliicker coordinates of boundary conditions 
of the problem (3) from its eigenvalues, if its differential equation is known. 

2. The Uniqueness of the Solution of the Inverse Problem. 

Together with problem (3), let us consider problem (4). 

Theorem 1 (on the uniqueness of the solution of the inverse problem). Sup- 
pose the following conditions are satisfied: 



If non-zero eigenvalues of problems (3) and (4) are identical, with account taken 
for their multiplicities, then problems (3) and (4) are also identical. 

Proof. A general solution of the differential equation of the problem (3) is the 
function y(r) = y(r, A) = C\ Jo(Xr) + C2Yo(Ar), where standard notations for 
cylinder functions are used (see [2]). 

The boundary conditions are used to determine the constants C\ and C<i- 
The following function is a characteristic determinant of the problem (3): 



(S). 



rank^4 = rank^4 = 2. 



(5) 
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A (A) ee det{AB T ), where B is a matrix of order 2x4 such that b\\ = Jo (A a) , 

&12 = Jo(Aa),6i 3 = Jo(A6),6i4 = Jo(Aa),6 2 i = Y (Xa) } b 2 2 = Y (Xa) } b 2S = 
Yo(Xb) , 624 = ^o(A6) . The equation for frequencies, is obtained from the condi- 
tion of the existence of a non-zero solution for C%.. The latter solution exists if 
and only if, the determinant A (A) is equal to zero (see [5]). 
Using the Binet-Cauchy formula, we obtain 

A(A)= A r Ai(A) = 0. (6) 

l<i<j<4 

It follows from asymptotic estimations for cylindrical functions 



J v [z) ~ \ — cos(z — utt/2 — 7r/4), Y v {z) ~ y — sin(z — utt/2 — 7r/4), 

\ 7T Z ' ' \ 7T Z 

that A (A) is an entire function (see [4]). Hence, it follows from Hadamard's 
factorization theorem (see [4]), that characteristic determinant A (A) of problem 
(3) and characteristic determinant A (A) of problem (4), are connected by the 
relation 

A(A) = C X k A(A), (7) 

where k is a certain non-negative integer and C is a certain non-zero constant. 
From this, we obtain the identity 

1<«<7<4 l<«<j<4 

It follows from this, that 

(Aij ■ Bij(X) - A l3 ■ C X k Bij(X)) ee 0. (8) 

l<i<j<4 

Note, that number k in this identity is equal to zero. Actually, let us assume the 
opposite: k 7^ 0. Using Maple, we get that the functions Bij(X) (1 < i < j < 4) 
and also the same functions multiplied by X k , are linearly independent. 

From this and identity (8), we obtain Aij = = 0(l<i<j<4), which 
contradict the condition (5) of the theorem. 

Hence, k = and 

J2 (Atj - C Aij) • ByW = 0. (9) 

l<i<j<4 

From this, by virtue of the linearly independence of the corresponding func- 
tions, we obtain A^ = C Aij } i.e. Pliicker coordinates of boundary conditions of 
the eigenvalue problems (3) and (4) are identical. 
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Then, it follows from Proposition 3, that problems (3) and (4) are also identical. 
The theorem is proved. 

3. Exact Solution of The Inverse Problem and Its Stability. 

This section is concerned with solving this problem and constructing an exact 
solution of the inverse problem. 

Suppose Ai, A2, A3, are the eigenvalues of problem (3). We substitute the 
values Aj, (i = 1, 2, 3) into (6) . Since Ai 2 = A34 = 0, we obtain a system of three 
homogeneous algebraic equations in the four unknowns ^13, A U} A 2 ^ } A 24l : 

F-Z T = 0, (10) 

where F = {fij\i = 1 , 2 ; j = 1,2,3,4} is a matrix of order 3x4 such that 
fn = B 13 (Xi), f u = B u (Xi), /i 3 = #23(Ai), «24 = B 2 4(\i), (i = 1,2,3), and 
Z = (A13, A u , A 2 s, A24)) is a row-vector. 

The resulting system has an infinite set of solutions. It follows from the unique- 
ness theorem, which has been proved that the unknown minors ^13, Au, A 2 % } ^24 
can be found apart from a constant. Hence, the resulting system must have a rank 
of 3 and a solution, determined apart from a constant multiplier. 

If the minors (Pliicker coordinates) ^13, Au } A 2 % } A 2 ^ are found apart from 
a constant, then then Problem 2 is solved. We can reconstruct any one of the 
equivalent matrix of boundary conditions by its Pliicker coordinates (see [3]). 

For example, if ^13 = 1, A\ 2 = A^ = then 

, 1 A 23 

1^! 

This reasoning proves 

Theorem 2. If the matrix F of system (10) has a rank of 3, then the solution 
of the inverse problem of the reconstruction boundary conditions (3), is unique. 

Note, that theorem 2 is stronger than theorem 1. Theorem 2 used only three 
natural frequencies for the reconstruction of the boundary conditions and not all 
natural frequencies as theorem 1, are used. 

It is significant that we have the following theorem on stability of the solution: 

Theorem 3. Suppose that one of the third-order minors of matrix F is sub- 
stantially non-zero. If \Xi — \i\ < 5 « 1, then the boundary conditions of 
problem (3), are close to the boundary conditions of problem (4)- 

The proof is analogous to the proof of the theorem on stability in [1]. 

4. Approximate Solution. 

Since small errors are possible when measuring natural frequencies, the problem 
arises to find an algorithm for the approximate determination of the type of 
fastening from the natural frequencies, which are found with a certain error. 



(11) 
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Let ^i4> ^23) ^24 be solution of the system of homogeneous algebraic 
equations(lO). It is unnecessary for these liquids to be totally immiscible. 

Aside from measurement and calculation errors, it is unnecessary for the values 
Al 3} A° lA , A 2S} ^4 to be minors of a matrix. So this problem is not trivial. We 
must find minors A13, Au } A 2 3, ^24, closer to the values ^3, A° 1A , A 23} A 2A . 

It is known from algebraic geometry [3] that the numbers A\2, ^-13, A\±, ^23, 
^24, ^34 are minors of some 2x4 matrix, if and only if, the following condition 
is satisfied: Ai 2 A M — Ai^A 2A + A^A^ = 0. This condition is called the Pliicker 
relation. 

If the Plucker relation for these numbers is realized, then Al 3 , A\ A , A^i ^24 
are minors of some matrix and corresponding boundary conditions are found. 

If the Plucker relation for numbers A° l2 , ^13, A° lA , A 23 , A 24} A^ 4 is not realized, 
the required minors A\ 2 , A13, Au, ^23, ^24, ^34 are found with the help of 
orthogonal projection. By definition, put x\ = A12, X2 = ^34, £3 = ^13, £4 = 
—A24, £5 = A14, xq = A23. Using this definition, we get the Plucker relation 

x\ X2 + X3 X4 + £5 xq = 0, (12) 

which characterizes a surface S in the 6-dimensional space. By definition, put 
I/i = A° 12 , y 2 = A° M , y 3 = A° 13 , y A = -A° 2A: y h = A° u , y 6 = A° 23 . 
By definition, put 

X = (x h x 2 , x 3 , X4, x 5 , x 6 ), Y = (y u 2/2, 2/3, 2/4, 2/5, 2/e), 
X* = (x 2 , x h x 4 , x 3 , ar 6 , x 5 ), Y* = (y 2 , 2/1, 2/4, 2/3, 2/6, 2/s), 

(X, Y) = Xi 2/1 + £ 2 2/2 + £3 2/3 + X 4 2/4 + ^5 2/5 + ^6 2/6, 



where X = OX 7 Y = OY , O is the origin. 

— * 

Let X be an orthogonal projection of Y on surface (12). The vector X* is 
normal for surface (12) at the point X. It is identical to the following equations 

Y = X + pX\ (13) 

where p is a real number. Having solved a set of linear equations (13), with the 
unknowns Xi, x 2 , £3, £4, £5, Xq, we obtain 

X = - ] —(Y-pY*). (14) 
1 — p z 

From (15), it is easy to obtain (see [1]) 



(Y,Y)-^{Y,YY-{Y,Y*Y 
p= - — • (15) 

(y,y*) 
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5. Example. 

If Ai = 2.93, A2 = 6.16, A3 = 9.34 corresponding to the first three natu- 
ral frequencies uji determined, using instruments for measuring the natural fre- 
quencies, then the solution of system (10), apart from a constant, has the form 
A° 13 = 1.00 C, A° u = 2.00C, A° 2A = -2.00C, A° 23 = -C., Using (14) and (15), 
we get A13 = 0.72 C, A u = 1.17 C, A 23 = -1.17 C, A 2A = -1.89 C. Suppose 
C = I/-A13; then from (11), we obtain 

1 -1.62 
"OOl 1.62 ' 

Note that the numbers Ai = 2.93, A2 = 6.16, A3 = 9.34 presented above are 
almost the same as the first three exact values, that correspond to a case with 

1-200 
12' 

This means that the unknown fixing of the membrane, which is inaccessible for 
direct observation, has been satisfactorily determined. The errors of determina- 
tion of boundary conditions, are stipulated with round-off errors of calculations. 
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